clc;
clear;
tic,
disp('------------------方法1-------------------------------')
f=@(x,y,z)z.^2.*exp(-(x+y.^2));

yM=@(x)sqrt(1-x.^2);
zm=@(x,y)sqrt(x.^2+y.^2);

zM=@(x,y)sqrt(2-x.^2-y.^2);

I=integral3(f,0,1,0,yM,zm,zM,'RelTol',1e-20);
toc
I

disp('------------------方法2-------超级慢，不建议运行------------------------')
syms x y z,
zm=sqrt(x^2+y^2);
zM=sqrt(2-x^2-y^2);
tic
I=int(int(int(z^2*exp(-(x+y^2)),z,zm,zM),y,0,sqrt(1-x^2)),x,0,1)
vpa(I)
toc
%历时 57.785078 秒。